Viscous Spacetime Fluid and Higher Curvature Gravity 
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The Einstein field equation as an equation of state of a thermodynamical system of spacetime is 
reconsidered in the present Letter. We argue that a consistent interpretation leads us to identify 
scalar curvature and cosmological constant terms representing the bulk viscosity of the spacetime 
fluid. Since Einstein equation itself corresponds to a near-equilibrium state in this interpretation 
invoking f(R) gravity for nonequilibrium thermodynamics is not required. A logically consistent 
generalization to include the effect of so called 'tidal forces' due to the Riemann curvature is pre- 
sented. A new equation of state for higher curvature gravity is derived and its physical interpretation 
is discussed. 
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Higher curvature Lagrangians have been discussed as 
logically possible constructions for gravitation since long 
1] , however current interest in such theories is inspired 
by various shades of quantum gravity/effective field the- 
ory 0- Amongst them f(R) gravity, where / is a non- 
linear function of the scalar curvature i?, is considered as 
an attractive model to explain the observed cosmic accel- 
eration phase 3. The conflict of the metric formulation 
of f(R) gravity with the solar system observations has 
been claimed to get resolved in the Palatini formalism 
[H, however it has been argued that this theory results 
in the appreciable deviations from the microphysics (e.g. 
electron-electron scattering experiments), and violation 
of the equivalence principle [f| • It may be asked: Is there 
an alternative approach to higher curvature gravity? 

Departing from the variational principle approach, re- 
cently a thermodynamical derivation of the f(R) gravity 
has been proposed |6[ generalizing the previous deriva- 
tion of the Einstein field equation Q • Thermodynamics 
of spacetime in [tJ is motivated by the black hole ther- 
modynamics: it is assumed that similar to the black hole 
entropy formula a universal entropy density a per unit 
horizon area for all local Rindler horizons could be de- 
fined, and the Clausius relation TdS = 5Q of equilibrium 
thermodynamics holds for all local horizons. Here the 
heat flow across the horizon SQ is defined to be the boost 
energy of the matter and T is the Unruh temperature. 
Jacobson then argues that since entropy is interpreted 
as horizon area the Clausius relation would be satisfied 
provided the spacetime curvature in the presence of mat- 
ter is such that the Einstein field equation holds. In Q 
this approach is applied to the assumed entropy density 
proportional to f(R), and it is found that local energy 
conservation entails entropy production term, and hence 
nonequilibrium thermodynamics of spacetime. Further 
the higher curvature equation of state obtained is shown 
to be identical with the field equation derived from the 
f(R) Lagrangian. 

Obviously rederivation of the known field equations by 
itself is not of much significance unless new insights are 
gained. Unfortunately the paradigm of field theory has 
overshadowed this promising thermodynamic approach 



in [1, 0|- Moreover we have shown Q that covariant 
divergence law for the matter stress tensor is satisfied 
in a conformally related spacetime, g^ v = f{R)gfiv with 
no need to postulate bulk viscosity production term as 
is done by Eling et al in [6|, and therefore the issue of 
nonequilibrium thermodynamics becomes trivial. Note 
that this result is independent of the unimodular per- 
spective suggested in Q. 

In this Letter spacetime as a thermodynamical system 
in the spirit of Jacobson's approach is considered with 
the two- fold aim: 1) to gain new insights into the na- 
ture of spacetime assumed to be some kind of fluid at 
large scale, and 2) to motivate nonequilibrium thermo- 
dynamics in the case considered in [7[ and generalize it 
to include the effect of the so called 'tidal force' caused 
by the Riemann curvature of spacetime @ and derive a 
new higher curvature equation of state. 

First we argue that a consistent equation of state inter- 
pretation of the Einstein field equation demands a careful 
reanalysis of the local equilibrium condition: the <E> term 
in Eq.(6) below while easily interpreted in the field the- 
oretic setting, seems to require a different meaning since 
local energy conservation is assumed to be a form of first 
law of thermodynamics [6|. This term in the Einstein 
field equation in the absence of matter stress tensor is 
suggested to represent a dissipative contribution to the 
stress tensor of vacuum. Analogous to the bulk viscosity 
for the dissipative fluid [lOj, <& term, therefore signifies 
viscous spacetime, and the Einstein equation already cor- 
responds to a quasi-equilibrium system. 

Next we recall Bondi's argument [ll| noted in [§[ that 
the essence of observable gravitation is that relative ac- 
celeration varies at each spacetime point. It would mean 
that even the assumption of local thermodynamic equilib- 
rium for local Rindler horizons would, in principle, fail. 
However the gravitation is so weak that the correction 
due to Riemann curvature is proposed to be incorporated 
retaining key elements of the Jacobson's approach in the 
following. The notion of local Rindler horizon is used to 
define heat flow and entropy change in Q. The equiv- 
alence principle enables locally flat spacetime at each 
spacetime point, and for an infinitesimal 2-surface ele- 
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ment through the point one can assume vanishing shear 
and expansion to the past (or inside) of the plane. The 
past horizon, the local Rindler horizon is assumed instan- 
taneously stationary. This idea helps in introducing an 
approximate local boost Killing vector field which can be 
related with the horizon tangent vector and afiinc param- 
eter along the null geodesic. Now instead of employing 
the Jacobi deviation equation for relative acceleration, we 
assume that the imprint of the tidal force is carried by 
the congruence of the null geodesies via a correction term 
dependent on the Riemann curvature. Both nonuniform 
temperature and viscosity lead to dissipation; since we 
do not know how to define these quantities here instead 
of the entropy change we calculate change in the thermo- 
dynamic potential F = ST. Heat flow and change in the 
horizon area are calculated for local Rindler horizons to 
lowest order in the affine parameter as is done in [Tj. A 
correction term proportional to the Riemann curvature 
is introduced in the change SF. Equating 8Q and SF, 
and imposing covariant divergence law for matter stress 
tensor finally give the desired equation of state. The 
Letter concludes with a brief discussion on the physical 
interpretation and prospects of this equation. 

Briefly the main ideas of 0] are as follows. The black 
hole formula, namely the proportionality between en- 
tropy and the horizon area, is assumed to hold for all local 
Rindler horizons at each spacetime point of the manifold 
M. Causal horizon at a point p is specified by a space-like 
2-surface B, and the boundary of the past of B comprises 
of the congruences of null geodesies. Assuming vanishing 
shear and expansion at p the past horizon of B is called 
local Rindler horizon. The energy flux across the hori- 
zon is used for heat energy, and calculated in terms of 
the boost energy of matter: define an approximate boost 
Killing vector field x M future pointing on the causal hori- 
zon, and related with the horizon tangent vector fc M and 
affine parameter A by x M = —a\k^. The heat flux to the 
past of B is given by 

SQ = J T^ifdZ" (1) 

The integral is taken over a small region of pencil of gen- 
erators of the inside past horizon terminating at p. If area 
element is dA then = k^dXdA, and Eq.(l) becomes 

5Q = -a J T^k v \d\dA (2) 

Change in the horizon area is given in terms of the ex- 
pansion of the congruence of null geodesies generating 
the horizon SA — J 6d\dA. The expansion of the null 
geodesies generating the horizon is given by the Ray- 
chaudhuri equation 

dO 2 

^ = — - - RvJPW (3) 

Assuming vanishing shear and neglecting 9 2 term we get 
the solution 

6 = -XR liV k> i k v (4) 



Assuming universal entropy density a per unit horizon 
area it is straightforward to calculate the entropy change 

SS = -a J R llv k^k v \d\dA (5) 

The condition that the Clausius relation is satisfied for 
all null vectors k^ , and making use of the Unruh tem- 
perature ha/2ir gives 

+ $g„» = {2v/na)T^ (6) 

The unknown function $ is determined using the covari- 
ant divergence law for the stress tensor and contracted 
Bianchi identity; Einstein equation with a cosmological 
constant (CC) A is obtained. Here Newton's gravita- 
tional constant is identified as G = l/Aha, and the func- 
tion <f> is obtained to be 

$ = ~ + A (7) 

In the light of thermodynamic derivation a logically 
consistent physical interpretation of the equation of state 
is proposed here. In the Clausius relation integrands of 
(2) and (5) have been equated for all null vectors to 
obtain Eq.(6). Since is the matter stress tensor, 
from the integral form itself it is logical to infer that 
R^u is proportional to the stress of assumed spacetime 
fluid. What does $ term represent? In the derivation $ 
term is needed to satisfy the first law of thermodynam- 
ics (in the form of local energy conservation) , and hence 
it should correspond to an additional stress tensor for a 
dissipative process in Eq.(6). Noting its formal similarity 
with the bulk viscosity stress tensor for an isotropic fluid 
proportional to the tensor <5y, see Eq. (8.4.42) in 10] it 
seems reasonable to identify $ term in Eq.(6) as bulk vis- 
cosity tensor of the spacetime. Deeper insight is gained 
analysing the CC term. Setting matter stress tensor zero 
and assuming a vanishing CC the Ricci flat spacetime 
becomes indistinguishable from G^ v = R^ v — \g^v = 
as R = 0. Thus vacuum Einstein equation with van- 
ishing CC could be interpreted as an equation of state 
of a perfect fluid. The presence of matter is suggested 
to cause a dissipative process leading to a nonvanishing 
scalar curvature, i.e. a bulk viscosity to the spacetime; 
geometrically scalar curvature arises via the contracted 
Bianchi identity so that the Einstein tensor satisfies the 
covariant divergence law. Here R has been given a ther- 
modynamic significance. 

What is the physical significance of A ? In the light of 
its appearance in combination with R in (7) this would 
correspond to the bulk viscosity treating the whole 4> 
term as the bulk viscosity stress tensor of the spacetime. 
Thus empty spacetime with nonvanishing CC would be 
like a viscous spacetime fluid. In Jacobson's approach the 
thermal behavior of quantum vacuum in flat spacetime is 
an important ingredient. The time translation symme- 
try in Rindler wedge corresponds to the boost symme- 
try of the Minkowskian spacetime; quantum vacuum in 
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Minkowskian spacetime as observed in the uniformly ac- 
celerating frame acquires thermal properties of a thermal 
bath of real particles with the Unruh temperature 
Assuming that quantum vacuum fluctuations would be 
significant only at Planck length scales a universal en- 
tropy per unit horizon area and Unruh temperature have 
been assumed in Q- This is a reasonable assumption, 
however in the light of thermodynamic derivation it may 
be asked: Does there exist averaged out observable ef- 
fect of vacuum fluctuations? There is an interesting phe- 
nomenon in quantum optics [l3| : atomic motion in light 
radiation in the so called optical molasses bears resem- 
blance with that of a particle in a viscous fluid. Analo- 
gous to this, envisaging gravitational molasses having ori- 
gin in quantum fluctuations in the spacetime bulk viscos- 
ity associated with CC could be attributed to this effect. 
Another possibility is, following suggestions in the liter- 
ature, to treat CC as the energy density of quantum vac- 
uum fluctuations, A = 8ttGVq/c where Vq is the vacuum 
expectation value. In this case, nonzero CC would cause 
dissipative process similar to matter energy and constant 
scalar curvature would determine the bulk viscosity of 
spacetime. We emphasize that both CC and quantum 
vacuum possess rather speculative character in general 
relativity, and the microscopic nature of the spacetime 
is unknown to us, hence the preceding discussion though 
quite plausible, also remains speculative. 

The question of bulk viscosity for the Einstein equa- 
tion is also discussed in [6(] ; the expected bulk viscosity of 
STia/Ait is obtained from the f(R) equation in the Ein- 
stein frame. Authors argue it to be incorrect, and note 
the ambiguity of sign. However negative bulk viscosity 
could be related with the acausal teleological boundary 
conditions for the black hole The stretched- horizon 
formalism developed in 14] gives useful insights on black 
hole physics; in particular, energy and momentum con- 
servation laws for the membrane resemble with those of 
viscous fluid. The present interpretation that the Ein- 
stein equation represents a viscous spacetime fluid could 
be viewed as a generalization, albeit a radical one, of the 
membrane paradigm of [14j . 

Thermodynamic derivation of the Einstein equation 
based on arbitrary spacelike 2-surfaces has been recently 
established [HI]. It is, therefore, reasonable to assert 
that the preceding physical interpretation is of general 
validity. Clearly viscosity term implies that the system 
is not in thermodynamic equilibrium; moreover higher 
curvature, e.g. f(R) gravity, is not needed to motivate 
nonequilibrium thermodynamics. The present interpre- 
tation however does offer an effective procedure to in- 
clude the effect of Riemann curvature as a correction 
to viscous stress tensor. Returning to the fluid anal- 
ogy, in general the dissipative coupling coefficient rjiju 
is a fourth rank viscosity tensor [101 ] . and has symmetry 
of indices similar to curvature tensor Rp Va p- Assum- 
ing that the correction due to the Riemann curvature is 
proportional to Rp aV p the product with the Ricci tensor 
keeping in mind the symmetry of the indices would be 



R ap R jlcJV p. Thus we construct the simplest form of 5F 

SF = -^1 ([R lll , + p{R a PR^ vp ~^R^)]k^k l '\d\dA 
In J 2 

The last term in the square brackett in the integrand 
corresponds to the viscosity as identified above in Eq.(7) 
for the scalar curvature. This also makes the application 
of covariant divergence law unambiguous. Using Clausius 
relation and equating the integrands of (8) and (2) for all 
null vectors we get 

Rpu + (i{R ap R^ vp - + *9w> = j^ T nv (9) 

The unknown function \l/ is determined taking the co- 
variant divergence of (9) and using the Bianchi identity 
and vanishing divergence of matter stress tensor. Use is 
made of following relations 

(JW^T = {\RopR ap -\R^),, + \{R;^ (10) 



(RR^ = {R^y + (4- - R-.a* ) ,f (11) 



The expression for W is finally obtained to be 



R I3R 2 PR^PR, 
if = 1- t 1 

2 8 4 



ap 



(12) 



Substituting ^ in (9) we get the desired higher curvature 
equation of state 

r _^CR - — n \j-RH<rP(Tl Ra Pn \-—T 

' ' (13) 

Perusal of the higher curvature field equations in the lit- 
erature [2|, y, [4[ shows that Eq.(13) is a new result. This 
equation has some remarkable properties: 1) The arbi- 
trary constant (3 could be adjusted to keep intact the es- 
tablished physical tests of the Einstein equation since the 
higher curvature effect appears as a correction to G pv . 2) 
The trace of (13) leads to the same result as that of pure 
general relativity. Note also that the first two terms are 
same as the ones in R 2 gravity in the Palatini formal- 
ism which would make it to be in agreement with some 
of the results of that theory. However the third term is 
a new addition. And 3) In contrast to the higher cur- 
vature field equations obtained from the action principle 
which arc fourth order derivative equations, Eq.(13) is a 
second order derivative equation similar to that one ob- 
tains in the Palatini version where the action is varied 
taking metric and affine connection as independent vari- 
ables. These characteristics strongly suggest viability of 
Eq.(13) as higher curvature gravity theory. 

It is well known that the question of gravitational en- 
ergy in general relativity is quite complex 16], and even 
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more so in higher curvature gravities [2|, llTj . An interest- 
ing consequence of (13) is that in the empty spacetime it 
reduces to 

+ pR°P(R^ vp - ^g^) = (14) 

What is the physical significance of the second term in 
Eq.(14)? In analogy with the traceless stress tensor of 
electromagnetic field it is plausible to identify this term 
to represent stress tensor of gravity radiation loss. Bondi 
presented a lucid discussion [l8[ on the inductive and 
wave transfer of gravitational energy in general relativ- 
ity, and considered axially-symmetric vacuum solution of 
Weyl and Levi-Civita for this purpose. Since the vacuum 
defined by Eq.(14) has different nature, it would be inter- 
esting to delineate inductive and wave energy transfers 
in this case. 



The problem of entropy production for f(R) gravity 
has been extensively discussed in [p|. What would be 
the entropy production in the present case? There is 
another interesting question: Could one derive Eq.(13) in 
the action formalism? We propose to investigate Palatini 
formalism for the action used by Deser and Tekin 
elsewhere since the main theme of the present work is 
thermodynamic approach to the spacetime. 

In conclusion, we have argued that Einstein equation 
in thermodynamic approach represents a viscous space- 
time fluid, and derived a new equation of state for higher 
curvature corrections. 
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